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Abstract 

In this paper we investigate the stability of a recently introduced Bose-Einstein con¬ 
densate (BEC) which involves logarithmic interaction between atoms. The Gaussian 
variational approach is employed to derive equations of motion for condensate widths 
in the presence of a harmonic trap. Then we derive the analytical solutions for these 
equations and hnd them to be in good agrement with numerical data. By analyzing 
deeply the frequencies of collective oscillations, and the mean-square radius, we hnd 
that the system is always stable for both negative and positive week logarithmic cou¬ 
pling. However, for strong interaction the situation is quite different: our condensate 
collapses for positive coupling and oscillates with hxed frequency for negative one. 
These special results remain the most characteristic features of the logarithmic BEC 
compared to that involving two-body and three-body interactions. 


Keywords Logarithmic BEC; Variational Method; Numerical Simulations; Linear 
Stability; Collective Oscillations. 

1 Introduction 

Since its experimental observation in trapped alkali atomic vapors m, Bose-Einstein 
condensation (BEC) has largely attracted much attention from both theoretical and ex¬ 
perimental viewpoints 111 [ 3 ] m 121. Among all the recent developments, the study of col¬ 
lective excitations occupies an important position in studying the properties of trapped 
BEC. Due to the diluteness of BEC, most theoretical works on collective excitations 
are mainly focused on considering the two-body interaction by means of the Gross- 
Pi taevskii equation (GPE) |l6l|7][8l|9l. A long side with the experimental progress with 
BECs in atomic waveguides and on the surface of atomic chips, which involve a strong 
compression of the traps, a signihcant increase of the density of BECs can be achieved; 
thus, three-body interaction may also play an important role in the process of collective 
excitations ifTOlfTTlfT^flBl . and stability of trapped BEC lfT4l [151 fl^ . Specially, it is 
reported in Ref. m that the addition of three-body interaction can extend consider¬ 
ably the region of stability for a condensate, and allow the number of condensed atoms 
to increase even when the strength of the three-body force is very small. However, this 
correction remains insufficient given the fact that in higher densities regime nothing 
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prevents the emergence of multi-body interactions. In this case, one should consider 
all the powers of 4’ being the condensate wave-function, which implies the usage 
of the non-polynomial functions that will necessarily appear in the wave equations for 
condensate. 

Recently, K.G. Zloshchastiev El introduced a new quantum Bose liquid as a 
candidate structure of physical vacuum, this structure is described by a non-linear 
Schrodinger equation of a non-polynomial kind, namely the logarithmic one ifTSl [T^ 
M- Taking into account small vacuum fluctuations, the author demonstrated by a sim¬ 
ple model that the generated masses of the otherwise massless particles can be naturally 
expressed in terms of both elementary electrical charge and extensive length parameter 
of the non-linearity. Later in Ref. ED, it has been proved that the condensation gov¬ 
erned by the logarithmic Schrodinger equation (LogSE) has several differences from 
that described by the GPE, particularly, it possesses the self-sustainability property: the 
logarithmic BEC tends always to form a Gaussian-type droplet- even in the absence of 
an external trapping potential. It has also been proved that the nature of elementary 
excitations depend strongly on the background density which changes the topologi¬ 
cal structure of their momentum space. In this paper, we attempt to explore more the 
properties of the logarithmic BEC, namely those related to collective excitations in a 
harmonic trap. Our work is facilitated by variational approaches. Using the solutions in 
the Gaussian form, the LogSE is transformed into a set of differential equations about 
some parameters that characterize the condensate wave-function. It is confirmed that 
the analytical solutions of condensate widths equations conform to numerical data. By 
performing a stability analysis, based on two requirements for instability which gives 
us the critical points for collapse and which states that the BEC becomes unstable if the 
frequency of collective oscillations is zero or alternatively the mean-square radius of 
the condensate wave-function tends to zero in finite time, it is shown that our conden¬ 
sate is always stable, and that the collapse may occur only for large positive logarithmic 
interaction parameter. Einally, a comparison between the behavior of the logarithmic 
BEC and the classical one with two-body and three-body interactions indicates that the 
former is more stable than the latter. 

The paper is organized as follows. In Sec.II, we present an analysis of the spher¬ 
ically symmetric case, derive the governing equation for the condensate width, and 
discuss stability criteria. In Sec.Ill, the same procedure is applied for the axially sym¬ 
metric case, which reveals the effect of anisotropy. In Sec.IV, we report the numerical 
results of the present investigation. In Sec.V, we highlight the system of two-body and 
three-body interactions. Einal section is our conclusion. 

2 Analysis of the Spherically Symmetric Case 

2.1 Variational Approach and Governing Equations 

To begin, let us consider a boson gas with a fixed mean number of particles N, mov¬ 
ing in an isotropic harmonic trap. If the particle density is considerably high and the 
temperature of the condensate is low enough, the dynamics of the Bose-Einstein con¬ 
densed atoms can be described by the logarithmic Schrodinger equation (LogSE)El 
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where ijj = ip{r, t) being the condensate wave-function which is normalized to N, m 
is the atom’s mass. The parameter b measures the strength of the non-linear interaction 
and a is needed to make the argument of the logarithm dimensionless. 

Given a non-linear Schrodinger equation like O it is quite natural to tray a power 
expansion of the non-linear part ln(a^|'!/'P), one will find the following 


ln(a3|V^|2) 





( 2 ) 


The lowest order in gives the familiar Gross-Pitaevskii equation which de¬ 

scribes only two-body interactions. By including the next term of the power one gets 
the so-called Gross-Pitaevskii-Ginzburg equation where also three-body interactions 
have to be taken into account (see for instance Ref. It is clear that the above 

expansion makes sense only when the density is near the special value 1 /a^ but in gen¬ 
eral it is not always the case, for this reason the logarithmic non-linearity arises as an 
alternative for describing the multi-body interactions. 

Now V = y (r) is a static harmonic oscillator with spherical symmetry given by 
V = uj is the radial frequency of the isotropic trap. By using dimensionless 

variables, r ^ lr,t I = ^Jhjrru^, we redefine the condensate wave-function 


as 


(3) 

such that 

poo 

47r / dr = Na^ = n , 

Jo 

(4) 


where n is the reduced number of particles for the system. Thus, the dimensionless 
equation corresponding to Eq. O can be rewritten as 








(5) 


where A is the dimensionless logarithmic parameter which describes the proportional 
relationship between the logarithmic coupling and the level spacing of the harmonic 
oscillator 

A = ^ . (6) 

nw 

The chemical potential ^ is given by the eigenvalue solutions of Eq. ©, with 
p0= [-V2+r2-2Aln(|^|2)]0. (7) 


In order to analyze the dynamics of logarithmic condensate in an harmonic trap, it 
is convenient to follow the variational ar)proachir71 l23]| . Furthermore, we shall seek the 
solutions of Eq. (|5]) in the Gaussian form ll24ll 


^2 

<f{r, t) = A{t) exp[- -b iP{t)r^ + ia{t)], 


(8) 
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where A{t) is the amplitude, ^{t) is the width, a{t) is the linear phase of the con¬ 
densate, and /3(f) is the chirp parameter. The reduced number of particles n is given 
by 

n = = constant. (9) 

The Lagrangian density generating Eq. Q in spherical coordinates is 

CM = - A|H»(1„(|^|>) - 1). (10) 

By using Eqs. ([8]) and ( fTOb . the effective Lagrangian density can be obtained. 


Leff = 47r 


r'^C{ip)dr 


= 713/2^2^3 ^ (5 _ 4in(A)) 
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where the overhead dot represents time derivative. 

Analyzing the corresponding Euler-Lagrange equations. 


dLeff d dLeff 


= 0 , 


dq dt dq 

where q stands for A{t), ^{t), a{t) and /3(f) respectively, we obtain 


( 11 ) 
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Eliminating d between Eqs. (fT4li and (fTSl l one obtains 


(14) 

(15) 

(16) 


1 2A , 

2/?e = ^-e-y-4/32^, (17) 

By combining Eqs. (fTbl l and (fTTl i we get the following second-order differential equa¬ 
tion for the evolution of the width 


^3 C 

Eq. (ITSl l illustrates the motion of a particle along the positive ^ direction in the effective 
potential E(.f) and can also be rewritten as : 


dV{0 

d^ 


(19) 
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Figure 1: The effective potential versus the width ^ for different values of loga¬ 
rithmic parameter A 


where 

^(^) = y + ^ + 2Aln(e), (20) 

The interpretation of Eq. (fTSl l is straightforward. The first term of the right-hand side of 
this equation corresponds to the attractive effect of the harmonic potentiel, the second 
is proportional to and is related to the dispersive effect caused by kinetic energy, 
and the third one is proportional to and coming from the non-linear logarithmic in¬ 
teraction. In order to understand the condensate dynamics, we must focus on the study 
of Eq. (fTSl l. At the first sight we can see that it is a non-linear differential equation. 
Therefore, though this equation tends to simplify the original problem of Eq. (Ull, it 
is still complex and non-integrable considering that the only conserved quantity is the 
Hamiltonian. A formal solution of Eq. (flSl) is obtained by taking the Hamiltonian H 
of the point particle 

after integrating we get 


t = 



_ ^ 


( 22 ) 


However, this solution doesn’t provide sufficient information about the condensate dy¬ 
namics. Hence we infer that the present problem can be solved only by considering 
some approximations which will be discussed below in more detail. 

Now let us analyze the nature of the effective potential of Eq. ( l20l i. Eig. 1 depicts 
the potential energy curves as a function of ^ for different values of the dimensionless 
logarithmic parameter A. We notice that the potential has a simple shape and contains 
for each case one local minimum, which means that the only reasonable motion of 
the system is periodic (anharmonic) oscillation around the minimum of the effective 
potential V{^). 
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2.2 Oscillations and Stability of the Logarithmic BEC 

Obviously Eq. (fTSl l shows that the internal dynamic represented by the evolution of the 
width of the condensate is affected by the logarithmic interaction between atoms. The 
equilibrium width which corresponds to the stationary states of the condensates can 
be calculated by setting the gradient of the potential Eq. (l20l) equal to zero 


dViO 


- 0: 


which yields 

1 2A 

Eq. (l24l can be rewritten in polynomial form as 


(23) 


(24) 


^^ + 2ACo'-l = 0, (25) 

Note that for given A only positive real solutions from Eq. ( l25l l represents physically 
realized equilibrium width, in our case the relevant solution is given by 


Co 





(26) 


Expanding Eq. (fTSl l around the equilibrium width, we can obtain the dynamical equa¬ 
tion of the width 

^ = fo + ^sin(wrf + 0). (27) 

where A, 9 are real constant and is the frequency of the collective oscillations (in 
units of uj) which is given by 



Erom Eq. ( |2^ one gets 

r 3 2A 11/2 

‘^r = \l + - 775 -- ■ ( 29 ) 

I- [(A2 + 1)^/2 - A]2 (A2 + 1)1/2 

Note that the solution Eq. ( l26b is stable only if the frequencies of collective modes are 
real, otherwise the solution is unstable. The frequency of collective oscillations, is 
among the relevant quantities usually used in literature to analyze the stability of the 
condensate. One of the requirements for instability of the BEC is that is zero ll24l . 
this condition gives us the critical points at which the system collapses. Erom Eq. (l29l l 
one can derive that these critical points are given byQ 

Ac = ±i. (30) 

^Note that for the BEC with two-body and three-body interaction the corresponding parameter of inter¬ 
action can be complex. The imaginary part describes the effect of inelastic collisions on the dynamics of 
BEC GllEbl. Nevertheless this interpretation doesn’t correspond with our case as it leads to non-physical 
complex values of the width. 


6 










a 

U)r 


b 

^0 


T. 

2 ( 

i; 

i( 



A 


A 


-10 


-5 


5 


10 


-10 


-5 


5 


10 


Figure 2: (a) The oscillation frequency ujr and (b) The stable equilibrium width 
versus the logarithmic parameter A 

However, this complex values don’t have any physical meaning in view of the fact 
that A is the dimensionless logarithmic parameter, which describes the proportional 
relationship between the logarithmic coupling and the level spacing of the harmonic 
oscillator that must take real values. Thus, our condensate remains always stable and 
the collapse will not take place in the case of small strength interaction. Fig. 2(a) shows 
the frequency of collective oscillations as functions of the logarithmic parameter. We 
notice that the frequency of collective oscillations presents a relative dependance on 
the logarithmic parameter A. In fact, we see that this frequency decreases with the 
decreasing of A and tends almost to 1.4 for a large negative value of A. 

The mean-square radius is another issue of stability, briefly speaking, collapse oc¬ 
curs as the mean-square radius of the condensate wave-function, (r^) = n~^ f r'^\ip\‘^dr, 
tends to zero in finite time, (see, e.g., EH EH), in the present units the corresponding 
expression is given by 



(31) 


taking into account that </> governed by Eq. Q realizes an extremum of the Hamiltonian 
and using the following scaling 


r 


4> = i ; with 


(32) 


where I plays the role of a Lagrange multiplier, it is easy to express the Hamiltonian 
and the chemical potential in term of (r^), for the Hamiltonian we obtain 



(33) 


where 




(34) 




( 35 ) 
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The minima of H are given by the roots of the identity 


5H{1) 

SI 


lf=i=0, 


then we obtain the following characteristic relation 


X, - (r2) - 3AZ, = 0 . 


( 36 ) 


(37) 


and the Hamiltonian becomes 

H^lXs + Ur^)-XYs. 

b D 

In the same way, we can obtain the corresponding chemical potential as 

= 2{r^} - X{2Ys - 3Zs). 


(38) 


(39) 


from Eq. OTl i one gets 

t^ = 3e- A(2y, - 3Z«). (40) 

From Eq. OTT i it is clear that an examination of the effect of the mean-square radius 
on the stability of the logarithmic EEC system, can be attained by considering the 
equilibrium width. Fig. 2(b) shows the equilibrium width as functions of the logarith¬ 
mic parameter A. We see that the equilibrium width decreases with increasing A and 
tends to zero for a large positive value, this means that it is only for a sufficiently large 
positive logarithmic interaction between atoms that collapse can take place. 


3 Axially Symmetric Case 


Analogous to the procedure of last section, we consider the logarithmic equation in 
dimensionless variables for an axially symmetric case. The governed equation for the 
logarithmic condensate wave function (p{p, z, t) at radial position p, axial position z 
and time t can be written as 


-.3 


19 / 9 ^ 

2p dp \ dp J 


1 9 ^ 1 

2^^ ^ 2 




-I- A?z^ I — A 2 In 




(41) 


where V{p, z) = p'^+X^z'^ is the axial trap. Here length and time are expressed in units 
ofl{= and uj~^ respectively, with w the radial trap frequency, Ai = lOz/uj being 

constant describing the anisotropy of the trap which denotes the ratio of the frequency 
along the 2 : direction loz to the radial frequency w, and A 2 = -^ is the dimensionless 
logarithmic parameter. 

We therefore try the product ansatz 


ip{p,z,t) = A{t) P[ exp[- ^ -f ide{t)g^ + 


2 


(42) 
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Figure 3; (a) The equilibrium width (b) The frequency Wz, versus the dimension¬ 
less logarithmic parameter A 2 for several values of Ai 


with normalization 


/ OO poo 

dz / p\Lp{p,z,t)\^ dp = Na^ = n , 

-00 J 0 


(43) 


The Lagrangian density in this case is 




(44) 


whereas the effective Lagrangian is given by 


/ OO poc 

dz / p£{(p)dp 

-OO J 0 


— + —pr + -— (5 — 41n(z4)) + -(dp + d^) 


' 4^2 ■ 2 ' 2^ " 


(45) 




(2/3p + API + 1 ) + I 


{Wz 


+ APi + \i 


As in the isotropic case one can write the Euler-Lagrange equations for ai{t), A{t), 
Pi{t), and After some algebra one can eliminate the variables ai{t), A{t), and 

Pi{t) from these equations and obtain the following second-order differential equations 
for the evolution of the condensate widths 


F -l_r 


■>P 

2 X 2 

e. 17 

Now the equilibrium widths equations are given by 




^PO “ t3 


1 

?p0 


2 A 2 

^pO 


(46) 

(47) 


(48) 
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(49) 


Al^zO 

and their solutions are respectively 


1 




zO 


2 X 2 

^zO 


^P0= (Al + l)^/'-A; 


-, 1/2 


(50) 


,5.) 

Using Eqs. (fSOl l and (fSTl i. we find that the frequencies of oscillations around these 
equilibrium widths are given respectively by 


Ulp — + 


ojz — Ai 1 + 


2 A 9 


1/2 


[(A2 + 1)1/2 _ A2]2 + 


3Af 


2 A 2 


1/2 


(52) 


(53) 


[(A2+Ai)i/2-A2P [^x\ + Xlf^-X2^ 

From Eq. (l53l l. we find that the values of A 2 which lead to the axial collapse are given 
by 

A 2 = ±/Ai. (54) 


Which are again non-physical complex values. Figs. 3(a) and 3(b) demonstrate re¬ 
spectively the equilibrium width and the frequency of axial mode as functions of the 
logarithmic parameter for different values of the anisotropy parameter. We see that the 
width decreases and the frequency increases with increasing either the positive loga¬ 
rithmic parameter or the anisotropy. In addition, if the interaction parameter is positive 
and very small, the anisotropy makes a great effect on both width and frequency, but 
when the interaction parameter becomes larger, this effect disappears and the width as 
well as the frequency are almost determined by A 2 . Figs. 3(b) shows also that if the 
negative value of A 2 gets larger, the frequency reaches a relative constant value depend¬ 
ing on the parameter of anisotropy Ai. Furthermore, with much larger anisotropy, the 
equilibrium width tends to zero, which practically represents a collapsed condensate. 


4 Numerical Results 

The main results of our investigation are set out in graphic forms. The graphs pre¬ 
sented here show the evolution of the condensate width in spherically (Figs. 4(a)-4(f)) 
and axially (Figs. 5(a)-5(d)) symmetric traps when Eqs. (fTSl) and (l47l i are solved by 
using a numerical routine. This numerical routine was implemented through the use of 
Mathematica for some different values of the logarithmic parameter A (or A 2 ) and the 
anisotropy parameter Ai. 

It is clear that these numerical results are in a good agreement with the analytical 
ones. Indeed, for the spherically symmetric case (Figs. 4(a)-4(f)) the condensate width 
performs in all cases small oscillations about the equilibrium width, furthermore, by in¬ 
creasing A the period decreases (the frequency increases), whereas, the period is equal 
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Figure 4: Evolution of the condensate width in spherically symmetric case (a) = 

0.951308 and A = 0.1, (b) (o = 0.643 and A = 1, (c) = 0.223 and A = 10, (d) 
^0 = 1-051 and A = -0.1, (e) = 1.55377 and A = -1, and (f) = 4.477 and 
A = -10 
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Figure 5: Evolution of the condensate width in axially symmetric case (a) Ai =0.1 
and A 2 = 0.1, (b) Ai = 0.1 and A 2 = -0.1, (c) Ai = 0.1 and A 2 = 10 (d) Ai = 0.1 
and A 2 = —10 


almost to 5.5 (the frequency ~ 1.14 ) if A = —10. Another aspect of this agreement 
appears in the value of the frequency itself, for example the frequency drawn from 
Eq. (l28l l for A = 10 is 28.39 and from Eig. 4(c) we can find that it is equal to 28.45. 

The same agreement is also observed for the case of axially symmetric trap, namely 
if we compare on the one hand Eig. 4(a) with Eig. 5(a) and Eig. 4(d) with Eig. 5(b), and 
on the other hand Eig. 4(c) with Eig. 5(c) and Eig. 4(f) with Eig. 5(d), we find that 
the frequency increases for both positive and negative interaction parameter with the 
increasing of the anisotropy, but when the positive interaction parameter is much larger 
the anisotropy has no great effect on the frequency. 


5 Bose-Einstein Condensate with Two-body and Three- 
body Interactions 


The present section reviews mostly the results of earlier works concerning the EEC 
with two-body and three-body interactions ll24l |28l and gives a comparison between 
this system and the logarithmic one. We should start writing down the governed equa¬ 
tion in a dimensionless form 




+ v + X2\^\^ + \M' 


'■P-. 


(55) 


In Eq. (I 55 I) we are assuming dimensionless variables; the unit of length is I = ^Ti/vnuj', 
and the unit of time is 1 /tu. V = V{r, z) is a static trap potential, that we assume a 
harmonic oscillator with axially symmetry given by E = -f Afz^), Ai = 
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is the anisotropy ratio. A 2 and A 3 are respectively the parameters of the two-body and 
three-body interactions which are considered reals. In particular A 2 is proportional to 
the two-body scattering length a^c and given by A 2 = Anasc/fv^l'^■ 

Next we give the analogs of the above constructions. Using the same apparatus of 
the last section, we write down the equations for the widths 


■ 

= 4 - 


p 


+ 


p 


K 

Wp 

K 


(56) 


(57) 


Here, we have introduced the dimensionless two-body and three-body interaction strengths 
given respectively by 


X2N dAaiV^ 

(2^)3/2 ’ ^ - 9V37r3 ’ 

where N is the number of atoms. The equilibrium widths are given by 


(58) 


1 . P 

c3 t 

?r 0 SrOSzO 


c3 '^ISzO + ^2 
SzO SrOSzO 


K 

jE~pr 

SrOSzO 


= 0 . 


K 


c3 

SrOSzO 


= 0 . 


(59) 

(60) 


Now let us analyze the stability of the present system with the emphasis on the simple 
case of isotropic trap. In this case we have similarity between Eqs. (l56l l and (l57l i. 
consequently the width equation, effective potential, equilibrium width equation and 
the frequency of collective oscillations are given respectively by 




1 P K 

^ ^ ^ 

Sr Sr Sr 


( 61 ) 




p 




+ 72 ) + 


K 


1 PK 

jr~^ro + jr + 7r-^- 

>r 0 ^rO ^rO 


(62) 

(63) 



where Eq. ( | 6 ^ was used in Eq. ( l64t . Therefore, for the repulsive two-body interaction 
case, we have only one solution for Eq. (l 6 Tb , then the condensate is always stable. As 
for the case of an attractive two-body interaction, we have a different situation: the 
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Figure 6: Condensate width = 'C^o = ■Co for (a) K = 0, as a function of P\ (b) 
K = 0.007, as a function of P; (c) P = —0.5, as a function of K. Solid blue lines 
represent the stable solution, dotted black lines represent another stable solution, and 
dashed green lines represent unstable solution. 


equation can have no equilibrium position, or it could have up to three equilibrium 
solutions. The results are summarized in Fig. 6 as the variation of the equilibrium 
width in term of the interaction parameters obtained by solving numerically Eq. (l6Tb . 
Fig. 6(a) shows solutions for K = 0 as a function of P. As we see, the solid blue 
curve represents the stable solution, whereas, the dashed green one depicts the unstable 
solution. In addition the critical point is identified as the join of the two curves, its value 
is Pc = —0.535 which coincides exactly with the value from Ref. fT). Between the two 
curves we have the stable area (S) such as if the initial width lies inside (S) the collapse 
is avoided, otherwise the collapse occurs (region U) even if the interaction strength is 
still finite. 

Next, considering the combination of an attractive two-body interaction and a small 
repulsive three-body interaction, we get different results of the stability. The system 
can either have one or three solutions, as shown in Fig. 6(b). Despite being small, the 
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addition of a positive three-body interaction K leads to the existence of at least one 
stable solution for each value of P. Nevertheless, as we can see from Fig. 6(b), for 
large negative values of P the curve representing ^o{P) is flat, which means that the 
stability region can be considerably enhanced by the inclusion of a small positive value 
of K compared to the case of pure two-body interaction. Furthermore, if the negative 
values of P is large enough, the equilibrium width tends to zero, which practically 
represents a collapsed condensate. 

In Fig. 6(c) we see that the stable solution for P = —0.5 exists only for a limited 
interval of negative values of K which leads to the fact that both attractive three-body 
and two-body interactions cause the same effect of instability on the EEC. 

Now let us explain the differences between the behavior of the EEC with logarith¬ 
mic interaction and the EEC with two-body and three-body interactions. Eirst of all, 
we begin with the most interesting cases of positive A and negative P (with K = 0 
), it is clear that both Eq. (fTST l and Eq. (|6TJ express otherwise the dynamics of one¬ 
dimensional soliton confined in harmonic trap. It is well known that in addition to the 
attractive term caused by the trap, such soliton objects are formed by a compensation 
between a dispersive and non-linear effects. In the two equations the dispersive effect 
is repulsive and proportional to but the non-linear effects are both attractive and 
proportional to for Eq. (fTST l and for Eq. (IbTT i. The remarkable difference be¬ 
tween the two systems will be more apparent if ^ 0. In fact, at this limit and for 

a small A the dispersive term dominates the non-linear one A^“^ for the soliton 
with logarithmic nonlinearity (Eq. (fTSl l). Therefore, as mentioned above the system 
remains stable and the collapse can take place only if A gets larger, i.e the non-linear 
term becomes larger. However, for Eq. (IbTT i the non-linear term P^~'^ dominates the 
dispersive one and the above compensation is no longer possible. Consequently, 
the EEC with two-body interaction collapses at the center of the potential. Yet, the 
addition of a small repulsive three-body interaction (depending on to the pure 

two-body one tends to push the soliton away from the center and to extend the region 
of stability, but even the stability region of the EEC with an attractive two-body in¬ 
teraction can be considerably extended by the addition of a small repulsive three-body 
interaction, it never reaches that related to the EEC with logarithmic interaction, and 
this fact remains one of the most discriminating features between the two systems. 

6 Conclusion 

The Gaussian variational approach is used to determine the dynamics of an hypothet¬ 
ical EEC governed by the so-called logarithmic Schrodinger equation. Despite being 
non-linear, the evolution equation for the width has provided a great simplification of 
the original problem, and its analytical solution has been verified to be in a good agree¬ 
ment with numerical results. This equation has demonstrated that the condensate width 
evolves like a classical particle in anharmonic potential. 

The stability of a trapped logarithmic condensate studied in the present work by 
considering two important quantities; the collective oscillations frequency and the 
mean-square radius. Eased on the relevant requirements for instability which gives 
us the critical points for collapse and which states that the EEC becomes unstable if 
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the oscillation frequency is zero or alternatively the mean-square radius of the conden¬ 
sate wave-function tends to zero in a finite time, it has been shown that our condensate 
in isotropic harmonic trap is always stable, and that the collapse is possible if the loga¬ 
rithmic strength is positive and sufficiently large. For the case of axially harmonic trap, 
it is proved that once the anisotropy gets larger, the equilibrium width tends to zero for 
either positive or negative values of logarithmic parameter. Consequently, beside the 
effect of logarithmic strength, there is also the effect of anisotropy which can promote 
collapse, in other words this phenomenon occurs also as the condensate takes a very 
thin pancake-shape. 

Finally, the comparison between the EEC with logarithmic interaction and the EEC 
with two-body and three-body interactions has definitely demonstrated that the former 
is more stable than the latter. 
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